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One possible channel for black hole formation is the collapse of a rigidly rotating massive neutron 
star as it loses its angular momentum or gains excessive mass through accretion. It was proposed 
that part of the neutron star may form a debris disk around the black hole. Such short-lived massive 
disks could be the sources of powerful jets emitting cosmological gamma-ray bursts. Whether the 
collapse creates a disk depends on the equation of state of the neutron star. We survey a wide 
range of equations of states allowed by observations and find that disk formation is unfeasible. We 
conclude that this channel of black hole formation is incapable of producing powerful jets, and 
discuss implications for models of gamma-ray bursts. 


A canonical mass of neutron stars born in supernova 
explosions is M ss 1.4M0. The distribution of M around 
1.4M0 might, however, extend above ^Mq, especially if 
the neutron star is born spinning fast, with a period ap¬ 
proaching the minimum (“breakup”) Pgpin 1 ms. The 
additional centrifugal support allows a stable hydrostatic 
configuration with mass M that would be forbidden for 
non-rotating stars. 

Neutron stars in binary systems have additional 
chances to gain mass through accretion. The second most 
massive known pulsar J1614-2230 is in a binary system 
and has M « 2 M 0 [1]. Its spin period is 3.15 ms. It is 
unclear if accretion is capable of spinning up the star to 
Pspin ~ 1 ms. If this happens, the star mass may keep 
growing and remain stable even when it exceeds 2Mq. 

Such centrifugally supported “supramassive” neutron 
stars (SMNS) may also be created in mergers of neu¬ 
tron star binaries. Recent observations of J1614-2230 
and J0348-I-0432 [1, 2] indicate that the equation of state 
(EOS) of dense nuclear matter is relatively stiff, and 
therefore some mergers may initially result in a stable 
object supported by pressure and fast rotation [e.g., 3]. 
Numerical simulations show that the object will initially 
rotate differentially [4-8], but that solid body rotation 
will be rapidly established following outwards transport 
of angular momentum via magnetic stresses and gravi¬ 
tational waves. The timescale for differential rotation to 
be removed could be as short as tens of ms [9], and will 
almost certainly be much shorter than 10 s (e.g., [10]). 
The heat stored in the merger product is also mostly lost 
to neutrino emission within seconds (e.g., [11]). 

The SMNS is fated to collapse to a black hole. Its 
lifetime is controlled by the eventual loss of angular mo¬ 
mentum (spindown-induced collapse) or excessive mass 
growth (accretion-induced collapse). The collapse is as¬ 
sociated with a huge release of gravitational energy and 
could produce a bright transient event — a burst of elec¬ 
tromagnetic radiation, such as a cosmological gamma-ray 
burst (GRB). 

This GRB trigger is plausible if the equatorial part 
of the neutron star is not immediately swallowed by the 
black hole but forms a compact, massive, centrifugally 


supported disk around it. Jets of hot plasma and radi¬ 
ation are expected to emerge from the debris disk and 
power the burst (e.g., [12]). 

In the merger scenario, the SMNS eventually collapses 
due to its gradual spindown, which removes the ro¬ 
tational support in minutes to hours. The spindown 
timescale depends on the magnetic field of the merger 
product, which is likely amplified to i? ~ 10^® G dur¬ 
ing the merger [13, 14]. This implies a moderate delay 
of the collapse-powered burst following the gravitational 
waves that are emitted during the merger and hopefully 
detected by Advanced LIGO [15, 16]. 

The goal of this Letter is to assess if the key condition 
for this burst scenario — a massive debris disk after the 
collapse — can be satisfied. The structure of the SMNS 
and hence the outcome of its collapse are controlled by 
the EOS of the dense nuclear matter P{p). Available gen¬ 
eral relativistic simulations of the collapse do not show 
disk formation [17, 18]. These simulations, however, im¬ 
plemented only simplified EOS. In particular, [17] used 
the polytropic P oc with index n < 2, and found 

that less than 10~^Mq remains outside the black hole at 
the termination of the simulation, comparable to their 
numerical resolution. They also found that for an ex¬ 
tremely soft EOS (with n = 2.9 and 3) disks can form, 
however such EOS are incompatible with observations of 
neutron stars. The remaining open possibility is that a 
different form of the EOS could lead to disk formation, 
e.g. soft at high densities (which gives a compact inner 
core — the seed for a future black hole) and stiff at lower 
densities (which gives an extended outer core with a high 
angular momentum). 

In this Letter we explore a wide range of EOS in search 
for one that could possibly give a debris disk. Instead of 
carrying out full-fledged and computationally expensive 
hydrodynamic simulations of SMNS collapse, we employ 
a simple method. We analyze the equilibrium hydrostatic 
configuration prior to the collapse and check if it satisfies 
a necessary condition for formation of a debris disk after 
the collapse. 

Condition for Disk Formation.— A stringent criterion 
on disk formation can be derived by assuming that all but 
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an infinitesimal amount of the SMNS’s mass and angular 
momentum are inherited by the newly formed Kerr black 
hole. Matter at the SMNS equator has the largest spe¬ 
cific angular momentum, jei and hence is the most likely 
to comprise a disk. The angular momentum is conserved 
during collapse, as long as magnetic and viscous torques 
are negligible and the spacetime remains axisymmetric. 
The centrifugal barrier will stop the equatorial matter 
from plunging the horizon if je exceeds the specihc an¬ 
gular momentum of the inner-most stable circular orbit 
(ISCO) in the Kerr metric of the nascent black hole, 

je > Jisco(a) ^ disk formation is possible. (1) 

Note that jisco depends on the spin parameter a = 
Jc/GM^ where J is the angular momentum inherited 
by the black hole from the SMNS. A similar criterion has 
been employed previously to the collapse of supermassive 
gas clouds [19]. 

Maximally Rotating Maximal Mass.— We construct 
axisymmetric neutron star models using the rns code 
[20, 21], which calculates relativistic rotating hydrostatic 
equilibria following the method outlined in [22, 23]. The 
collapse occurs when the stellar mass exceeds M^ax at 
which the star becomes secularly unstable [24] and no 
hydrostatic solution is found. 

Afniax depends on the angular momentum J and the 
EOS of dense nuclear matter. For a given EOS, we cal¬ 
culate Mniax(</) and find a and je immediately prior to 
collapse. Disk formation is clearly impossible for a non¬ 
rotating star because matter will fall radially into the 
newly-formed Schwarzschild black hole. As J and hence 
je are increased, black hole spin a increases and hence 
jisco decreases. Condition (1) could thus in principle be 
satisfied at some point along the maximal mass sequence. 

The maximal mass sequence cannot be ex¬ 

tended indefinitely as it eventually reaches the mass- 
shedding limit, beyond which the co-rotating orbital fre¬ 
quency at the SMNS equator exceeds the SMNS rota¬ 
tion frequency. This point defines the maximally ro¬ 
tating maximum mass (MRMM), Mniax(Anax), which is 
typically 10-30% higher than M„iax(0). The collapsing 
MRMM has the best chance to form a debris disk but 
this is not guaranteed. Although je of the MRMM is 
just sufficient to orbit the hydrostatic star, the space- 
time metric changes after the collapse and the same je 
can fail to sustain Keplerian rotation around the nascent 
black hole. If condition (1) is not met for the MRMM, 
it will not be met for any slower rotating maximal mass 
models and we may conclude that disk formation is im¬ 
possible for this EOS. 

The input parameters of the rns code are the central 
energy density and the oblateness of the star. For a given 
oblateness we find the maximal mass model by varying 
the central energy density. Then we step along the max¬ 
imum mass sequence toward MRMM by increasing the 



FIG. 1. Sequence of neutron star mass M and spin param¬ 
eter a for three sample EOSs, illustrating our method for 
assessing the possibility of disk formation following SMNS 
collapse. The top portion of the figure shows the maximal 
mass sequence (triangles/squares/circles) and mass shed lim¬ 
its (small points) for each EOS. Masses are normalized to the 
maximal value for a non-rotating star corresponding to each 
EOS. The bottom portion shows the (dimensionless) specific 
angular momentum of a test particle at the SMNS equator, 
jec/GM, along the maximal sequence curves. A solid red line 
denotes the minimal angular momentum required to orbit the 
resulting Kerr black hole with spin parameter a, jisco (a). Ac¬ 
cording to the criterion (1), disk formation is ruled out as 
long as je lies below this red curve. These three EOS are also 
marked in Figure 2. 

oblateness parameter. At the end of the sequence we it¬ 
erate the oblatness until the mass shed limit is found to 
within a specified accuracy. At each step we check if the 
disk formation criterion (I) is satisfied. 

Figure I illustrates our procedure for three represen¬ 
tative EOS, labeled eosA, eosB, and eosC, respectively. 
For eosA, je < jisco for all black hole spin a, so disk for¬ 
mation is impossible according to criterion (1). For eosC, 
je > jisco for a > 0.5, indicating that a disk could form; 
however, the maximum non-rotating mass for this unre¬ 
alistically soft EOS is only O.48M0. Disk formation is 
also possible for eosB, but only for a very narrow range 
of J near the mass-shedding limit. 

Survey of the EOS Space.— The possibility of disk for¬ 
mation is controlled by the high density EOS, which is 
poorly known. Therefore, below we conduct a survey over 
a broad range of EOS. Our goal is to check whether it is 
possible to simultaneously satisfy the disk formation cri¬ 
terion and current observational constraints on neutron 
star radii and masses. 

We parametrize the EOS a.t p > po = g cm“^ as 

a broken power law. This choice is motivated by previous 
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FIG. 2. Regions of allowed and forbidden disk formation in 
the EOS parameter space. Dashed purple curves show con¬ 
tours of constant maximum mass for non-rotating neutron 
stars, while dotted black lines indicate constant radius values 
for a 1.4Mq non-rotating star. The green region shows the 
2cr allowed parameter space based on observed neutron star 
masses [2] (bottom boundary) and constraints on neutron star 
radius [27] (left and right side boundaries). 

works [25] which show that a piecewise polytrope can 
reliably reproduce a variety of EOS models. The break 
is fixed at density pi = 10^'^ '^ g cm“^. At densities below 
po we use the SLy EOS [26] with the approximation of 
[25], and we fix P(po) to the SLy value. 

With fixed pi we are left with only two free param¬ 
eters: Pi = P{pi) and the power-law index at p > pi, 
r2 = dlnP/dlnp. Two degrees of freedom in the EOS 
may be insufficient to predict observables to within ~ 1% 
accuracy (e.g., as in [25]). However, this form of EOS 
is sufficiently flexible for our purposes, allowing indepen¬ 
dent variation of the SMNS mass M and radius R. These 
parameters determine the star’s compactness M/R, the 
key factor for disk formation. 

The results of our numerical survey of the parame¬ 
ter space Pi-r2 are shown in Eigure 2. For “stiff” EOS 
above the grey strip even the MRMM configuration fails 
to meet the criterion (1), and thus disk formation is ruled 
out. The criterion is met by the MRMM below the grey 
strip (and possibly inside the strip where it is numerically 
unresolved). 

Small Pi or r2 values are however problematic as they 
predict low M^ax while observations demonstrate the ex¬ 
istence of neutron stars with M « 2Mq [1, 2], even at 
moderate rotation when centrifugal effects may be ne¬ 
glected (the 39 ms spin period of J0348-I-043 is slow 
enough that it can be treated as essentially non-rotating 
for the purpose of constraining the maximal neutron star 
mass). 


An additional observationaly accessible parameter is 
the radius of normal neutron stars with moderate rota¬ 
tion and canonical mass M ss 1.4M0. For instance using 
observations of transiently accreting and bursting neu¬ 
tron stars [27] reported RiaMq = 10.42 — 12.89 km at 
2(7. We note that current neutron star radius constraints 
are subject to uncertainties in both astrophysics and nu¬ 
clear physics modeling and the radius constraints are not 
entirely settled yet (cf. e.g., [28, 29]). 

For any candidate EOS one should check its predic¬ 
tion for Mniax(>^ ~ 0) as well as i?i. 4 MQ, which can be 
tested against observations. Figure 2 shows the contours 
of constant = 0) and RiaMq on the Pi-r 2 plane 

together with the observational constraints. The condi¬ 
tion Mmax > 2M0 alone excludes almost the entire re¬ 
gion where disk formation is possible. A significant gap 
appears between this region and the allowed region if fol¬ 
lowing [27] we also require RiaMq < 13 km. We also 
note that even below the grey strip formation of a debris 
disk requires significant fine-tuning toward the MRMM 
configuration. Disk formation quickly becomes impossi¬ 
ble if J is reduced below Jmax (see Figure 1, in particular 
model eosB). 

Discussion and Astrophysical Implications.— Our 
method employs a simple parametrization for the high 
density EOS as a piecewise polytrope, and hence may 
not replicate nuances of realistic EOSs (for instance, [25] 
shows that it tends to overestimate the sound speed). 
This parametrization is, however, sufficient to capture 
the overall mass distribution of the star, which is most 
important to our analysis. Since our results in Figure 2 
show a significant gap between the observationally al¬ 
lowed and the disk forming regions, any analysis using 
a more complex parametrization (as in e.g. [25, 27]) is 
likely to yield a similar conclusion. This is especially 
so when considering that much of the formally allowed 
disk formation region will not produce a disk in practice 
without fine tuning of the SMNS angular momentum. 

In particular, we have applied our method directly to 
the entire list of parametrized EOSs given in [25] (their 
Table III), and find that none of these support disk for¬ 
mation. Interestingly, the robustness of our main con¬ 
clusion relies in part on the recent discovery of a 2Mq 
neutron star [1, 2] and hence could not have been made 
with as much confidence prior to 2010, when the largest 
known mass was 1.74±O.O4M0. 

Although lower limits on the maximum neutron star 
mass are well established by dynamical measurements, 
observational constraints on the neutron star radius are 
subject to systematic uncertainties [e.g., 30]. It is thus 
important to note that our conclusion that disk forma¬ 
tion is unlikely depends most sensitively on the estab¬ 
lished maximum mass constraints, and less critically on 
the neutron star radius. 

Our analysis assumed axisymmetric collapse. This 
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is reasonable since non-axisymmetric perturbations will 
likely be damped out via gravitational waves. Further¬ 
more, if the amount of surviving disk mass is determined 
by deviations from axisymmetry then producing a disk 
of an interesting mass > \Q~^Mq translates into a radial 
perturbation of > 2 km, an unlikely occurrence. 

We have additionally assumed that magnetic or vis¬ 
cous torques do not affect the SMNS matter during the 
collapse. Numerical hydrodynamical simulations consis¬ 
tently show that the SMNS matter collapses on a dynam¬ 
ical timescale with approximate conservation of angu¬ 
lar momentum and negligible dissipation effects on fluid 
streamlines [17]. Magnetic fields could become dynam¬ 
ically important only when they are extremely strong. 
Such fields could also slightly affect the SMNS structure. 
Its radius would be increased up to ~ 16% in the most 
extreme case of magnetic pressure equal to thermal pres¬ 
sure (e.g., [32]). 

Our results have implications for some GRB models. 
Electromagnetic emission from SMNS formed in neutron 
star binary mergers has been proposed by many authors 
[e.g. 33-36] to explain long-lived X-ray flares (“extended 
emission”) and plateaus observed following short dura¬ 
tion GRBs [e.g. 37-39], which in some cases have been 
observed to terminate abruptly in a way suggesting a 
SMNS that has collapsed to a black hole [40]. These 
magnetar models have been criticized because it is not 
clear how to produce the relativistic jet responsible for 
the initial GRB itself as the result of baryonic pollution 
from the young neutron star remnant [e.g., 41]. This 
has recently led to the suggestion of a “Time Reversal” 
scenario [15, 16], whereby black hole formation and the 
GRB is delayed for tens or hundreds of seconds follow¬ 
ing the merger, but due to light time travel effects is 
observed before X-rays from the SMNS remnant cease. 
A similar physical situation, which posits the collapse of 
a SMNS to a black hole following the accretion of mat¬ 
ter from a binary companion (accretion-induced collapse; 
e.g., [42, 43]) is also commonly invoked as an alternative 
to neutron star merger models for short GRBs. 

Both these alluring models (accretion-induced collapse 
and Time Reversal) require a debris disk after the SMNS 
collapse in order to power the short GRB. Our results 
show that this assumption contradicts the stiff nuclear 
EOS inferred from observations of neutron stars. 

This does not necessarily mean that SMNS collapse 
will have no observational electromagnetic signature. For 
instance [44, 45] suggest that if the SMNS is initially 
magnetized, a significant electromagnetic transient could 
arise regardless of any surrounding accretion disk. How¬ 
ever, such a transient is unlikely to last many dynamical 
times across the black hole horizon and hence may fail to 
explain the 0.1 — 1 s duration of observed short GRBs. 

Our model assumes solid body rotation and a cold EOS 
and hence does not rule out a disk if the black hole forms 
shortly following a binary neutron star merger. Disk for¬ 


mation in fact appears to be a robust outcome of gen¬ 
eral relativistic simulations of the merger process [e.g., 
9]. Thermal pressure is only sustained for a few seconds 
after the merger, until neutrino cooling sets in. More 
importantly, the merger remnant is primarily supported 
by differential rotation, such that the collapse is usually 
initiated by the outwards redistribution of angular mo¬ 
mentum, as is expected to occur on a timescale of tens 
or hundreds of milliseconds due to magnetic or viscous 
stresses. Since in this case collapse occurs prior to the 
establishment of solid body rotation throughout the rem¬ 
nant, disk formation is much more likely than in the case 
of a delayed collapse. 

Finally, our results also render untenable proposed sce¬ 
narios for long duration GRBs which postulate a long de¬ 
lay (exceeding hours or days) between the core collapse 
of a massive star and the formation of a black hole with 
a debris accretion disk [46]. 
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